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The form of the vacuum stress-tensor for the quantized scalar field at a dielectric-vacuum 
interface is studied. The dielectric is modeled to have an index of refraction that varies with 
frequency. We find that the stress-tensor components, derived from the mode function 
. . . expansion of the Wightman function, are naturally regularized by the reflection and trans- 

, mission coefficients of the mode at the boundary. Additionally, the divergence of the vacuum 

■ energy associated with a perfectly reffecting mirror is found to disappear for the dielectric 
mirror at the expense of introducing a new energy density near the surface which has the 
opposite sign. Thus the weak energy condition is always violated in some region of the 

55 ' spacetime. For the dielectric mirror, the mean vacuum energy density per unit plate area in 

a constant time hypersurface is always found to be positive (or zero) and the averaged weak 
energy condition is proven to hold for all observers with non-zero velocity along the normal 
' direction to the boundary. Both results are found to be generic features of the vacuum 

(N ■ stress-tensor and not necessarily dependent of the frequency dependence of the dielectric. 
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5 ' I. INTRODUCTION 

o, 

, The renormalized vacuum stress-tensor for quantized fields at perfectly reflecting boundaries has been 

■ extensively studied, sufficiently so that they have made their way into various texts. Fulling and Birrell 
and Davies ||^ are two examples. The perfect reflector enforces a Dirichlet boundary condition that must 
be satisfied by the quantized field. For the case of the scalar or electromagnetic field with an infinite planar 

^ i' boundary, the renormalized vacuum energy density for stationary observers has the form 

(D ■ 

^ • (P)-. = ^, (1) 



k>( \ where a is a unitless constant of order unity, x is the transverse distance from the planar boundary and 

. TO > 1 is the dimension of the spacetime. The most notable property of the energy density is the divergence 

' in the limit a; — > 0. 

The divergence is not at all unexpected. The application of the Dirichlet boundary condition fixes the 
value of the field at a specific point with no uncertainty in the position. Thus the uncertainty of the conjugate 
momenta must be infinite, giving rise to the stress-tensor divergence. The precise mathematical origin of the 
divergence comes from the renormalization of the vacuum energy via the Green's function. The argument is 
simply that the unrenormalized Green's function is constructed such that it vanishes on the boundary. To 
renormalize the Green's function, the corresponding Minkowski space Green's function must be subtracted 
away. Since this is divergent in the limit a; — *■ 0, the subtraction carried out for renormalization will leave 
behind a divergent part. When the vacuum energy is calculated from the renormalized Green's function, it 
will diverge at x = also. 

The divergence of the vacuum energy density causes considerable consternation since the vacuum would 
contain an infinite amount of energy. In a self consistent theory of gravity, the back reaction of the vacuum 
energy on the spacetime would be considerable. Yet we know from experience that metallic foils do not cause 
significant gravitational distortions around themselves. 
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There are three resolutions to the problems: (1) Replace the Dirichlet boundary condition with something 
more realistic, which is the topic of this paper. (2) Add position uncertainty to the mirror; this has been 
discussed by Ford and Svaiter (3) Treat the entire problem quantum mechanically, which has been 
partially touched upon by Barton and Eberlein |^ and more recently by Heifer 1^. 

Replacing the Dirichlet boundary condition has to some degree been discussed by other authors. For 
example Candelas Q has derived the asymptotic form of the vacuum energy density for the problem of two 
nondispersive dielectric half-spaces which are in contact with one another along a planar interface. More 
recently Heifer and Lang Q carried out similar calculations for a nondispersive dielectric to vacuum interface, 
with considerable interest placed upon the energy conditions in the vacuum region. 

This paper addresses the behavior of the stress-tensor for a quantized massless scalar field at a planar 
dielectric to vacuum interface where the index of refraction of the material varies as a function of frequency. 
The Dirichlet boundary condition is replaced with appropriate frequency dependent reflection and trans- 
mission coefficients. The stress-tensor in the vacuum region is well behaved, even on the interface, and is 
largely independent of the functional dependence of the index of refraction. Therefore the vacuum energy 
divergence is a result of the Dirichlet boundary condition and not a generic feature of the field theory. It 
should be noted, however, that it is possible to get divergences for dielectric materials at the interface. The 
simplest example is a material with a constant index of refraction. In this case, the divergence comes about 
because part of the amplitude of the wave at every frequency still satisfies Dirichlet boundary conditions. 
The existence of the divergence is intimately tied to the failure of specific integrals to exist. For example, in 
two dimensions, a divergence will occur if 



fails to exist. Here R{lli) = [/^{uj) — 1] / [n{uj) + 1] is the reflection coefficient for the modes at the interface, 
while n{uj) is the index of refraction. There are similar integrals in four dimensions involving both the 
reflection and transmission coefficients. 

There are two remarkable results regarding frequency dependent dielectric models. First, all of the integrals 
that must be evaluated to find the components of the stress-tensor are self regularized, meaning that they 
all have a built in cutoff function that prevents ultraviolet divergences. For example, in two dimensions the 
energy density for a static observer is 



The only constraint placed on the index of refraction is for it to tend to unity as the frequency grows. Thus 
the reflection coefficient plays the role of the regulator in the integral. Similar results are found for the four- 
dimensional case as well, where there exists another integral for the evanescent modes which is regularized 
by the transmission coefficient. 

The other unique property, which we shall prove in general, is that the vacuum stress-tensor satisfies the 
averaged weak energy condition (AWEC) for any moving geodesic observer, provided that the component 
of the observer's four velocity normal to the dielectric to vacuum interface is non-zero. This occurs even 
though the pointwise weak energy condition (WEC) over some regions of the spacetime fails. The AWEC 
results from the Riemann-Lebesgue lemma applied to integrals of the form of Eq. (^) evaluated in the limit 
X ^ oo. In order to apply the lemma, it is required that conditions like exist. 

A related result is that the energy density, integrated along the normal direction in a constant time 
hypersurface in the vacuum region, is a positive constant. In two dimensions the constant is zero and 
independent of the choice of model for the index of refraction. Thus the spacetime always contains two 
regions, one where the vacuum energy is positive and another with equal negative energy. In four dimensions, 
the constant is found by evaluating the integral 



for a specific choice of n{uj). We consider only the class of models for which n{uj) > 1 for all frequencies. 
The function space for which the integral holds is obviously much larger, however this is beyond the focus 
of this paper. 




(2) 




(3) 




(4) 
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Specific models are examined in this paper. The vacuum stress-tensor is found to asymptotically approach 
the form of the perfectly reflecting vacuum stress-tensor in the large distance limit. The price of removing 
the divergence of the stress-tensor by using a frequency dependent dielectric is the introduction of a vacuum 
energy near the surface with opposite sign. This simultaneously is the reason why the AWEC holds while 
the WEC fails. 



II. MODE STRUCTURE AT A DIELECTRIC INTERFACE 



We begin with the source free scalar wave equation in m-l-1 dimensional Minkowski spacetime for a 
nondispersive material, 

(V2-n252^(/)(x,i) =0, (5) 

where n is a constant. Units of ?i = c = 1 are being used throughout. A positive frequency plane wave 
solution to this equation in an infinite medium takes the form 

/k(x,i)cxe'("'^-''-"*), (6) 

where lo = |k|. The modes can be normalized by defining the probability density for the wave equation and 
setting its integral over all space to unity. This leads to the normalization condition 

(/k(x,i),/k-(x,i)) = -* j d"a;nM/k(x,i)[a,/^,(x,t)]-[ajk(x,i)]/^,(x,t)}-^™(k-k'). (7) 

Thus the normalized plane wave mode takes the form 

/k(x, t) = n("-2)/2 [2cj(27r)"]"^/^ ^^(n\.■^-u>t) (g) 

The most general solution to the wave equation (j^) is then given as a mode function expansion, 

/ + 00 
d^k ak/k(x,t)+aj,/^(x,i) . (9) 



If 0(x, t) is a classical solution to the wave equation, then ak and are the Fourier coefficients for the mode 
function expansion. The classical solution has the property that it can in general be written as 

(t){'K,t) = g{x~-t) + h{x+-t). (10) 
n n 

This represents two traveling waves, g and h, that propagate in the +x and —x directions, respectively. Each 
traveling wave maintains a constant spatial profile that is just translated along the direction of motion. When 
second quantization is applied the cocfhcients Ok and aj^. become the annihilation and creation operators, 
respectively. 

Next, consider wave propagation in a dispersive material. A scalar wave equation is needed which mimics 
the behavior of the electromagnetic field inside a dispersive material. In electromagnetism, the wave equation 
is found by first Fourier expanding the Maxwell equations for harmonic time dependence. Upon combining 
them, we find that each component of E and B must individually satisfy the Helmholtz wave equation 

[V^ - n^{uj)uj^] (/)(x, Lu) = 0. (11) 

We will take this to be our dispersive wave equation. The mode functions which satisfy this equation are 
the same as given above with the replacement n — > n(w). The general solution (/)(x, i) is given as a mode 
function expansion in terms of the plane waves in the same form as (|^). To obtain this form we must 
impose the condition n*{uj) = n{uj). However, even a simple resonance model (anomalous dispersion) for 
the dielectric as given in Jackson |^ leads to an index of refraction that does not satisfy this condition and 
would necessitate a more complicated quantum field theory. Thus, we restrict n(uj) to be real for simplicity. 



3 



The physical interpretation of this condition is that the dielectric does not attenuate the wave propagating 
through it. It also excludes materials that may have any degree of conductivity. Strictly speaking, causality 
and the Kramers-Kronig relations require any index of refraction to have a nonzero imaginary part. However, 
it is also possible to have broad frequency bands within which the real part is very large compared to the 
imaginary part. Thus, for the models studied here, we are assuming that the dominant contribution to the 
energy density come from regions of the spectrum where Re{n) ^ Im{n). The general solution to the wave 
equation in the time domain no longer has the physical property of traveling waves as in the nondispersive 
material. For this case, each frequency component propagates with a different velocity, so the overall profile 
of the wave changes as it moves. 

Now consider the example of two different dielectrics which have a planar interface at x\ = 0. To the 
left of the planar boundary the index of refraction will be denoted by nj. To the right of the boundary 
the index will be ut- In addition we will define the normal to the interface, u = (1,0,- • - ,0). We allow 
a monochromatic plane wave to be incident on the interface from the left with a momentum k/. At the 
interface this incident wave is partly reflected back into the first medium and partly transmitted into the 
second medium. The reflected and transmitted waves have momenta and k^ respectively, and the three 
components of the wave are 

(/)„c.(x,0 =e*("^'^^™*) (12) 

(t>re,fl. (x, t) = R e«("^ k„.x-a;Bt) (^3) 

0(,a„.(x,t) =re^(''--™*) (14) 

For electromagnetism, the plane wave modes at the boundary would have to satisfy the continuity of the 
transverse and perpendicular components of the fields across the interface. For the scalar field, the boundary 
conditions that will be applied are the continuity of <p{x,t) and dx(l){^,t) across the boundary. This leads to 
several relations; 

1. The continuity relations must hold at all times, therefore we must impose 

= UJR = U)T = w. (15) 



2. All of the terms must have the same functional dependence on the surface of the interface, therefore 

n/k/ • X = riijkfl • x = kx • x at a;i = 0. (16) 

Taking the first two of these relations yields the law of reflection, 9i = 6r. Taking the first and third 
of these relations yields Snell's law, nisinOj = ut sin 9t- This allows us to solve for the components 
of kii and ky in terms of the components of k/, 

ki = {ki,k2,- ■ ■ ,km), (17) 

kfl = (-A:i,A:2,---,fcm), (18) 

kT = {/3,nik2,--- ,nikm), (19) 

where 

= ^n^w^-n] {kl + --- + kl). (20) 

For n/ > Ut, it is evident that j3 will be imaginary for some set of incident momenta. These are the 
evanescent modes for which the incident wave undergoes total internal reflection. At the interface, all 
of the incident energy will be reflected while there exists an exponentially decaying tail in the medium 
that has a lower index of refraction. This theory of scalar dielectrics simulates all of the properties for 
wave reflection and transmission of the electromagnetic field. 
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3. The result from the continuity relations are the reflection coefficient, 



R = — \ (21) 















2ki 













and the transmission coefHcient 

T= — (22) 



These relations are identical to the reflection and transmission coefficients for the transverse magnetic 
(TM) modes of the electromagnetic field. 

The mode functions can now be easily written down. There are two types, "right-going" which have ki > 
and "left-going" which have ki < 0. The right-going modes have the form 

{m-2)/2 

/k(x,t) = ^ {e(-a;) re^"^'=i^i + i?^„e-'"^'=^^il + e(a;)T^„e''^^^} e'"-'('==^=+-+'='"^"')e-''"*, (23) 



where Q{x) is the step function 



In addition the right-going reflection and transmission coefficients are Rrg = R and Trg = T as defined 
above. 

The left going mode functions have a similar form, 

(m-2)/2 

[2w(27r)">]''^ 

(25) 

with left-going reflection and transmission coefBcients given by the interchange of nj with ut- In addition, 
the change in sign of the x-component of the momentum must be accounted for. Thus the left-going reflection 
coefflcient is 



Ri9 = \ \ (26) 





{kl + -- 






(fc| + .. 




-2fci 







and the transmission coefficient 

Tig = , (27) 

-ki + \j{^) uj^-{kl + --- + kl) 

In addition, the new propagation constant is 

/?' = ^JnjLJ^-n^^{k^ + ■■■ + kl). (28) 

The normalization of the mode functions has been chosen such that the incident plane wave portion of 
each mode satisfies the orthogonality relation in an infinite medium. We are now prepared to begin our 
calculation of the vacuum stress-tensor outside of a frequency dependent dielectric material. 
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III. TWO-DIMENSIONAL DIELECTRIC HALF-SPACE 



A. Vacuum stress-tensor 



We begin with a two-dimensional space with a frequency dependent index of refraction, n = n{u)), occu- 
pying half the spacetime and vacuum in the other half: 



n ■■ 



n(ui) for a; < 0, 
1 for a; > 0. 



(29) 



There arc two simplifications to the mode structure in two dimensions. Both result from the lack of 
spatial dimensions transverse to the interface. First there are no evanescent modes, thus the mode structure 
is characterized solely by fci = fc with u = \k\. The positive frequency right-going modes, k > 0, then take 
the form 



fk{x) 



1 



\/47ram 

Similarly, the left-going modes, fc < 0, are 



fk{x) = 



The second simplification is that there is just a single set of reflection and transmission coeSicients, 

R = and T - ^ 



n+1 " n + 1' 

respectively, which satisfy 

R^ + T'^ = 1. 

The Wightman Green's function for two points in the vacuum region of the spacetime is 

n+oo 



(30) 



(31) 



(32) 



(33) 



/-l-oo 
-oo 



1 

An 



" dk 



^ik(x-x') _ cos k{x + x') + i?2 e-*'=(^-^') 



-iu!(t-t') 



+ 00 



dk 



rp2 ^ik{x — x')^—iu}{t- 







1- 



In Minkowski spacetime, the Wightman function is 

'■+°° dk 



gife(a;-x')g-jai(t-t')_ 



(34) 



(35) 



We can renormalize the Wightman function in the vacuum half space by subtracting the Minkowski spacetime 
Wightman function to yield 



^Re„. {x, x') =D+{x,x')- i?+ {x, x'), 



1 



dk 



-2R cosk{x + x') + R^ e-*'^^^-^') 



-iu}{t—t') 



r+°° dk 
Jo ^ 



pa ^ik{x-x')^-iuj{t-t') 
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where we have made use of Eq. ( |33| ) . This can be further reduced by making the change of variable in the 
first integral of fc = — cj, resulting in 



1 



2tt Jq lo 

The stress-tensor in the vacuum region for arbitrary coupling is [^j 



(36) 



(37) 



where ^ is the coupling constant. The renormalized vacuum expectation value of the stress-tensor is found 
directly from the renormalized Wightman function 0, 



(0|rtt|0)He„. = lim 



(38) 



Special care is taken in ( pq ) to ensure that the derivative operators are symmetric under interchange of 
the primed and unprimed variables. While this does not affect the calculation of the energy density, it has 
serious effects for the remaining components of the stress-tensor. Substituting Eq. (^6|) yields 



(0|T„|0)n,„ = ^ / div LJ R{lj) cos2ljx. 

Jo 

Similar calculations show the remaining components of the stress-tensor vanish, 



1 




(39) 



(40) 



It is interesting that for ^ = 0, which is both the minimal and conformal case in two dimension, we have a 
vanishing stress-tensor. The form of the stress-tensor given in Eq. (^) requires this because of the fact that 
the trace must vanish for conformal coupling. That is, the only way the stress-tensor can have this form and 
have a vanishing trace is for it to be identically zero for ^ = 0. 



B. Energy Conditions 

Before we look at specific examples, we examine the general behavior of the energy conditions. First we 
look at the integrated energy density (total energy) in a constant time hypersurface. We define the energy 
contained within a volume bounded on the left by the dielectric surface and on the right by an imaginary 
boundary at a; = L as 

E{L)^ {0\Ttt\0)dx = -- R{uj)sm{2Luj)duj. (41) 
Jo I" Jo 

This yields the interesting result that in the limit L ^ oo the total integrated energy is zero. This follows 
directly from the Riemann-Lebesgue lemma for the ordinary Fourier integral The requirement is that 
/q R{uj)dio exist, which in general is true for models where we require limt^^oo ^(t^) = 0. Thus for the 
dielectric mirror, if the integral condition exists then there must be equal amounts of positive energy and 
negative energy. One can then infer that there is no vacuum energy divergence on the surface of the dielectric. 
It is interesting to note that the failure of the integral condition on the reflection coefficient corresponds to 
divergent vacuum energies, as demonstrated in the constant index of refraction example below. 

Not much can be said about the WEC other than it is necessarily violated at some time along a geodesic 
observer's worldline. Only by specifying the reflection coefficient can any specific comments be made. How- 
ever, it is remarkable that the AWEC can be proven in general. We begin with a geodesic observer whose 
world line is given by 
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x-(r) ^ (1 - ^ ^ u^{r) = (1 - v^'^' ^ ) , 

where r = [0,cxd). The vacuum energy density integrated along this worldhne is 



{T^^)ui'u''dT ^ --(l-v^)-^/^ hm / sin 2(1 - v^y/^tovr 

IT T^r- ' 



(42) 



(43) 



which, as we have seen above, is zero for any observer that has w > 0. Thus the AWEC is exactly satisfied. 
This is because such observers sweep through all of the energy distribution, encountering equal amounts of 
positive and negative energy along the way. For stationary observers, the AWEC may or may not be satisfied 
depending on the position at which the observers sit. 



C. Examples 

1. Constant Index, Vacuum and Perfect Reflector 

Let us consider the trivial case where n(uj) = constant = tiq for all frequencies. This makes the reflection 
coefficient, 

i?(^) = i?o = (44) 
no + 1 

independent of frequency as well. Thus the energy density is 

(0|rtt|0)H„„ = -Ro / dujLJC0s2ujx. (45) 

Jo 

However this integral is not well defined, and requires a regularization scheme to be evaluated. One method 
is to introduce a cutoff function, carry out the now well defined integral, and then take the limit as the cutoff 
goes to one for all frequencies. An example cutoff function is /(w) = e~°'^ in the limit as a ^ 0. Therefore, 
the regularization of the vacuum energy results in 



|r„|o)H„. = lim rdiu 



u!e 



TT [(2x)2 +a2]2' 



2Trx^ 



(46) 



The first thing to note is the divergence of the vacuum energy density at the surface of the dielectric. The 
divergence is expected from the Riemann-Lebesgue lemma, as the integral of the reflection coefficient over 
all frequencies is divergent. This is identical, up to the factor of Rq, to the known result for the perfectly 
reflecting mirror. Actually the index of refraction in this case is only modifying the numerical coefficient of the 
functional form associated with the perfectly reflecting mirror. Therefore, the replacement of the perfectly 
reflecting mirror with a dielectric results in the strength of the vacuum energy density being governed by the 
index of refraction. Low indices yield very weak vacuum energy densities relative to the perfect reflector, 
while high indices yield appreciable fractions of the perfect reflector vacuum energy. 

At one end of the spectrum for Eq. (^ ) is that of the dielectric having an index of refraction of unity. This 
is the case of the vacuum spacetime everywhere. In this case the reflection coefficient vanishes and there 
is no vacuum energy density. This result is consistent with that expected from vacuum renormalization in 
two-dimensional Minkowski spacetime. 

At the other end of the spectrum is the perfect reflector, which is deflned as the reflection coefficient being 
unity for all frequencies, R{i-u) — 1. Associated with this is an inflnite index of refraction which is rather 
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unphysical or at least ill defined. Thus, it may be better the think of the perfect reflector to be the limit of 
the reflection coefficient going to one, 

(0|rtt|0)p„fe,, Hcflcoto, = (47) 

This is the standard result for the vacuum energy in a two-dimensional spacetime with perfectly reflecting 
boundary conditions 



2. Discrete Cutoff 

The next example to consider is a dielectric which has a constant index of refraction up to a cutoff frequency 
LUc and unity thereafter, 

If we insert this into Eq. (|39|), we find 

(0|rtt|0)H.„. = -— f ^^^^ ) / dLuujcos2ujx, 
TT \no + lJ Jo 

where x = 2lOc x is a dimensionless quantity. From this example we see that the reflection coefficient is acting 
as the regulator to obtain finite results from what would otherwise be an ill defined energy density integral. 
Thus more realistic boundary conditions leads to a stress-tensor that is self-regularized. 

Eq. ( |49| ) is plotted in Fig. |l| where the most noticeable fact is the disappearance of the divergence in the 
vacuum energy density as the dielectric-vacuum interface is approached. The divergence has been replaced 
with a finite value of opposite sign. From the knowledge of the energy conditions above we know that this 
is a general feature of the vacuum energy for frequency dependent indices of refraction. 



3. Exponentially Decaying Reflection Coefficient 

An example of a frequency dependent index of refraction that does not have a discrete cutoff is 

n{Lu) = ^ . (50) 

We define a cutoff frequency ojc = ct~^ as the point where the reflection coefficient R{lOc) — 1/e. The 
associated energy density is 



(0|rtt|0)Hc„ = ^ / dujuje-°"^ cos2ujx, 

Jo 



2£ d_ 

TT da 

2^ CO,' X'-I 







2- (51) 

This function is plotted in Fig. ^ which clearly demonstrates the two characteristic behaviors of the vacuum 
energy for frequency dependent dielectrics. The first is that the vacuum energy asymptotically approaches 
the form of the perfectly reflecting vacuum energy at large distances from the mirror. The second is the 
introduction of a new energy density near the mirror surface which goes to a finite value on the surface and 
has the opposite sign to the perfectly reflecting vacuum energy. Integrating Eq. ( pT| ) over x = [0, 00) yields 
zero, confirming that there is equal positive and negative energy in the vacuum. 



9 



^. Gaussian Reflection Coefficient 



Finally, consider the case where the reflection coefficient is a Gaussian. The index of refraction is given by 



n{Lu) = 



2 2 



(52) 



As was the case in the preceding example, the cutoff frequency is defined as = a ^. The energy density 
is easily calculated. 



TT 



dujuje °' ^ cos2lox, 



l_7^|e~(*)'Erfi 



(I) 



(53) 



where Erfi(a;) = — iErf(ia;) is the imaginary error hmction. The behavior of the energy density is plotted in 
Figure |^, and we see that its behavior is quite similar to that for the exponential reflection coefficient of the 
preceding example. 



IV. FOUR-DIMENSIONAL DIELECTRIC HALF-SPACE 
A. Vacuum stress-tensor 

The dielectric to vacuum half space is now considered in four dimensions, with 

n{uj) for xi ~ X <Q, 
1 for xi — X > Q. 



(54) 



The mode structure has the form of Eqs. (^3|) and (25) and their respective transmission and reflection 
coefficients as derived in Section The renormalized Wightman function for both spacetime points in the 
vacuum region is 



1 



167r3 
1 

167r3 
1 

16^ 



dk-. 



— OO J —OQ 



° dk 



dkz I dk,, 



dk-. 



dk. 



-OO ^ 

° dh 

-OO ^ 

°° dk 



Rio 



j[ky{y-y')-\-k^(z-z'}-uj(t-t')] 



„|7;,^|2gi[(/32:-/3*a;')+nfc„(iy-y')+«fe.(2-^')-t^(t-*')] (55-) 



It is a straightforward but tedious task to act on the Wightman function with the appropriate derivative 
operator to find all the components of the stress-tensor. For simplicity we will consider the stress-tensor of 
the minimally coupled scalar field only. Three components are found to be non-zero. 



|T^.|0):,_ =diag[(0|rtt|0),0,(0|Ti,|0),(0|Tii|0)] 



(56) 



The formal expression for each is 



Ttt\0) = — ^/ uj^du) 2/ d6l sin^6'iT!i„ cos(2wa;cos6') -I- 



7r/2 
sin^ ^ {1/n) 



dO sh\''e\Trg\^ e 



2 -2ux^n^ sin^ 6-1 



(57) 
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and 



1 1 f°° \ 

(0|T;i|0) = --(0|Ttt|0) + — J / uj^duj 2/ dO sin 6 Rig cos {2ujx cos 

2 Stt-^ 7o 7^/2 



t/2 



+n 



7r/2 



at' smt^ \ lrn\ e ^ 



(58) 



where we have transformed to spherical momentum coordinates about the normal to the interface, 

{kx, ky, kz) {locos 9, wsinflcosc/j, a; sin 6* sin (p) , (59) 

and the integral over (p has already been carried out. The needed reflection and transmission coefhcients 
now take the form 



sm 



— COS! 



+ - sin^ 



and 



2 COS 6* 



cosy + %\ sm 



(60) 



The integral expressions involving the cosine and left-going reflection coefficient is the four-dimensional 
equivalent to that found previously for the two-dimensional case. The new term involving the right-going 
transmission coefficient is the contribution to the vacuum polarization due to the evanescent modes. Recall 
that evanescent modes do not exist in two dimensions, so we cannot relate these new terms to anything 
seen previously. However in higher dimensions, their importance to the vacuum stress-tensor cannot be 
understated. For example, in the energy density term, we see from the form of Eq. (57) that the evanescent 



modes contribute only positive energy density. The components of the stress-tensor are plotted as a function 
of position in Figure ^ for a specific choice of the index of refraction. 

Again we find that the stress-tensor components asymptotically approach the case of the perfect reflector 
in the large distance limit. Close to the dielectric to vacuum interface, their behavior is again modified taking 
on the opposite sign and going to a finite value. In the n oo limit for all ui, the dielectric stress-tensor 
reduces to the standard form of the perfect reflecting plate stress-tensor M 



|T^.|0)He 



1 



-diag[l, 0, -1, -1] 



(61) 



It should be noted that the perfectly reflecting plate stress-tensor is a Lorentz tensor that can be formed from 
the metric ?7p^ and the unit normal to the plate. However this is not a property of the dielectric stress-tensor 
in Eq. (^6|). In general, (0|Ttt|0) ^ — (0|r;;|0) for an arbitrary choice of the function n{uj). The failure of 
(0|Tp,y|0)Re„, to be a Lorentz tensor results from n{uj) being a frame dependent quantity. 



B. Energy Conditions 

Again, the mean energy density per unit plate area, the WEC and the AWEC are considered. We begin 
with the energy density as seen by a static observer integrated along the normal direction to the interface 
in a constant time hypersurface. We evaluate 



















J TT 



7r/2 



-i(l/n) 



sin^6l|T^g|2g-2c^:rV^ 



^6* 



(62) 



Here E represents the total amount of energy in a rectangular column, with unit base area at the dielectric 
to vacuum interface and of infinite length in the normal direction. Interchanging the orders of integration 
in the above expression, then evaluating in x gives 
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1 f"^ 



(63) 



where Ii, I2 and I3 are the three remaining angular integrals With the aid of Eq. (60), these integrals can 
be individually considered. The first is 



Ii = 2 Ri„ sin (2ll!Xq cos 9) dO 

Jn/2 cose 

M^O^^Oy) ^ 2^ r 2)^(^2 3in(2^^^y)rfy + 

y [n^ - 1) I Jo 



+ [2y{y^-l) + {n^-l)y]sm{2ujxoy)dyY (64) 

where the change of variable y = —cos0 has been made. At this point we apply the Riemann-Lebesgue 
lemma to the second and third terms of the expression when we take the limit xq — > 00. It is easy to show 
that 



{l-y^)^(n^-l)+y^\dy^^ 



-2n{n^ - 3) + (ji^ + 2n^ - 3) In ^ 



n + 1 
n- 1 



(65) 



and 



/ 2 1^1/2 1^1, / 5(n^ - + 5) for 1 < n < V3, , . 

\2y{y -l) + {n -l)y\dy^^^^ _ ^ forn > VS ^^^^ 

exist for all values of 1 < n < cx3. Therefore, both terms in Ii vanish, resulting in 

lim /i=-2 lim / sin(2a;a;oy) ^ _^ ^.^ S'i(2cjxo) = -tt, (67) 

xq — >oo xq — 'oo Jq y xq — ^00 

where Si{x) is the sine integral. 

The next angular integral to be evaluated, 

h = -n^ ^ /"'/^ |T;<,pe-W»^-n^e-id(?, (68) 



i-i(i/ri) V «^ sin^ 9-1 

has exponential suppression in the large a;o limit. The only place of concern would be when to — 0. However, 
the factor of tu^ in the frequency integral would cancel any contribution from that point. Therefore the 
entire integral vanishes in the limit xq 00. An exact analysis using Laplace's method for two-dimensional 
integrals to find the asymptotic expansion of this integral for large xq yields the same result. 
The final integral, 

h^-rv^r" , |r.,pdg = ^(n^ + 3), (69) 

Jsin-i(i/n) V« sin^ 9 -I * 

results in a constant that is independent of xq. When all the angular integral contributions are added 
together, the total energy contained in the column of unit area is 

/ Lj"^ \n^{uj)-l] duj. (70) 

Immediately, we see that the energy is always greater than zero for any model that has \ < n < 00. In 
addition, this term will be finite for models where n{ijj) 1 in the — > 00 limit. Therefore, such models 
will not have surface divergences in the stress-tensor. 

Since the vacuum stress-tensor asymptotically approaches that of the perfect reflector in the large distance 
limit, we can assume that the negative energy region is distant from the mirror, and most likely very small in 
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overall magnitude. Because of Eq. (|7(]|), the positive energy must be of significantly greater magnitude and 
close to the dielectric mirror. This behavior is evident in Fig. ^. However, until the frequency dependence 
for the index of refraction is given, nothing exact can be said about the WEC. 

The case of the AWEC is strikingly different. In fact, the averaged energy along the worldline of a moving 
observer with nonzero normal velocity is always positive definite. Consider the geodesic 



T 

VxT 

VyT 



where 7 = (l 



'1/2 



(71) 



and Vx ^ 0. We are assuming that the observer is moving outward from the mirror, so Vx is positive. The 
energy density averaged along the observer's worldline is 



{T^Ar))n..u^u''dT^-f^ lim 



[{0\Tu\0) + ivl + vl){0\Tu\0)]dr. 



(72) 



The first term in the integral is identical, up to factors of 7 and Vx, to that done to arrive at Eq. (^y). The 
same analysis using the Riemann-Lebesgue lemma on the second term involving (0|T;;|0) finds that nonzero 
contributions come only from the — (0|rtt|0) part when integrating Eq. (pq). Thus, we obtain the AWEC, 



1 



7 



1287r Vx 



UJ 



[n^iuj) - 1] dLj> 0. 



(73) 



The AWEC is satisfied for any outward or inward traveling geodesic observer's worldline. This results from 
the observer crossing the region of positive energy close to the dielectric to vacuum interface. However, 
for observers traveling parallel to the interface the AWEC can fail if the observer is in the negative energy 
region. Such geodesies never cross the positive energy region of the spacetime. 



V. SUMMARY 



In the preceding sections, we have shown that the general characteristics of the stress-tensor in the vacuum 
region outside of a planar dispersive dielectric can be discussed without specific knowledge of the frequency 
dependence of the dielectric. The divergence of the stress-tensor at a planar interface is linked to the failure of 
a set of integrals involving the reflection and transmission coefficients to exist. If such integrals do exist then 
the mode function expansion of the vacuum stress-tensor is self regularized. At sufficiently large distances 
from the dielectric mirror, the vacuum stress-tensor asymptotically approaches that of the perfect reflecting 
mirror. However, near the dielectric mirror the form of the stress-tensor approaches a finite value of opposite 
sign to that of the distant vacuum energy. This leads to the remarkable results that the mean energy per unit 
plate area is positive definite (or zero) and the AWEC is satisfied along all inward or outward half-infinite 
timelike geodesies. Taking the geodesic's four-velocity to unity yields the ANEC as well. 

The results here are for the free quantized scalar field. However the mode structure and boundary con- 
ditions are identical to that of the transverse magnetic modes of the electromagnetic field, so we would 
conjecture there would be comparable results for the electromagnetic field case. However the divergence 
that would hopefully be removed is not in the stress-tensor but in the expectation values of the square of 
the field strengths. 

There is a rich line of future topics to research with respect to the vacuum stress-tensor with more realistic 
boundary conditions. The most immediately interesting would be to carry out a similar analysis for the 
quantized electromagnetic field at a planar dispersive dielectric boundary. Part of the foundation for this 
has already been carried out by Heifer and Lang 0. In addition, a more realistic model for the index of 
refraction should be considered that includes conductivity and attenuation. The study of curved surfaces 
should also yield interesting results. One problem that stands out is the Casimir force (and stress-tensor) 
between two parallel infinite planar dispersive dielectric half spaces separated by a gap. This problem has 
an extremely interesting mode structure which includes evanescent, tunneling and waveguide modes. One 
can also consider mixing the dispersive dielectric model with non-zero position uncertainty for the dielectric 

i- 
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FIG. 1. Spatial variation of the vacuum energy density for the discrete cutoff in two dimensions. Unhke the 
perfectly reflecting boundary vacuum energy, given by the lighter curve, the energy density for the discrete cutoff 
does not diverge as one approaches the interface. Here the energy density is expressed in units of ( "°T? ) ■ 




FIG. 2. Spatial variation of the vacuum energy density for an exponentially decaying reflection coefficient in two 
dimensions. Near the dielectric mirror, the functional form of the vacuum energy density differs significantly from 
that of the perfectly reflecting mirror, given by the lighter curve. However, in the large x limit, the dielectric mirror 
vacuum energy density reproduces the behavior of the perfectly reflecting mirror. The energy density has been 
normalized to units of . 

7r 
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FIG. 4. Numerical simulation of the Tu and Tu components of the vacuum stress-tensor, solid and dashed lines 
respectively, in the free space region for a planar frequency dependent dielectric to vacuum interface in four di- 
mensions. The wave equation for the scalar field has minimal coupling and the model of the index of refraction is 
n(aj) = 1-1- 3 exp(— it)'^). The lighter line is the case of the perfect refiector. 
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